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Abstract. - We study a variant of the generalized Thomson problem in which n particles are 
confined to a neutral sphere and interacting by a l/r'' potential. It is found that for 7 < 1 the 
electrostatic repulsion expels all the charges to the surface of the sphere. However for 7 > 1 
and n > nc{-y) occupation of the bulk becomes energetically favorable. It is curious to note 
that the Coulomb law lies exactly on the interface between these two regimes. 



In a recent paper pP Bowick et al. studied a system of particles confined to the surface of 
a sphere and interacting by a repulsive 1/r'^ potential with < 7 < 2. They called this "the 
generalized Thomson problem" . It is interesting, however, to recall that the original Thomson 
problem was posed as a model of a classical atom Thus, n electrons were supposed to 
be confined in the interior of a sphere with a uniform neutralizing background, the so called 
"plum pudding" model of an atom. The Thomson problem, which is still unsolved, is then to 
find the ground state of electrons inside the sphere. 

In the absence of a neutralizing background the electrostatic repulsion between the particles 
"dynamically" drives the charges to the surface. This significantly simplifies the calculations 
by reducing the search of the ground state from the three dimensions down to two [3] . But 
what if instead of the Coulomb potential electrons interacted by a l/r"^ potential? Would 
they still go to the surface or prefer to stay in the bulk? This question was not addressed in 
the paper of Bowick et al. who have a priory confined their particles to reside on the surface. 

It is clear that for a small number of charges, mutual repulsion will force them to the 
surface. What happens, however, as the concentration of particles increases? To answer this 
question we compare the electrostatic energy of the configuration in which all n particles are 
on the surface of a sphere with a configuration in which n — 1 particles are at the surface and 
one particle is located at the center of a sphere. The electrostatic energy of n particles of 
charge q interacting through a generalized Coulomb potential g^/er^, with dielectric constant 
e, confined to the surface of a sphere with radius a can be obtained by considering the elec- 
trostatic energy of the two dimensional one component plasma (OCP) Ff^^'^ , i.e. charges on 
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the surface of a sphere with a neutrahzing background ^ISHSIE] , 

„2 9I-7 „2 01-7 

2ea'>' 2-7 ea'' 2 - 7 ^ ' 

The first term is the electrostatic energy of mutual repulsion between the charges, the 
second term is the self energy of the neutralizing background, and the third term is the energy 
of interaction between the charges and the background. The advantage of working with the 
one component plasma is that its ground state energy can be estimated by considering the 
interaction of an individual charge with the background inside its Wigner-Seitz cell. The 
characteristic distance d between the charges on the sphere is such that, ircpn = Ana? , and 

d-^. (2) 

The ground state energy of a Wigner crystal of charges interacting by a potential is, 
therefore, 

F^^^ = -M,n^ , (3) 

where is the Madelung constant. Of course, there is no perfect crystalline order of charges 
on the surface of the sphere and some topological defects must be present [S]. Nevertheless, 
we expect that the topological defects will not modify the scaling form of the Eq. ij^ll but will 
only affect the numerical value of the Madelung constant. Equation can now be rewritten 
as 

„2 r 91-7 _ Af . „, 1 

(4) 



En = ^ 



V?^ ' 



7 2-^-1 

This is equivalent to Eq. (4) of the reference T*. For the case of Coulomb interaction, 
7 = 1, it has been known for some time 9 that Eq. with Mi k, 1.102 gives an excellent 
approximation to the ground state energy of electrons on the surface of a sphere. This value of 
Ml is, indeed, very close to the Madelung constant of a planar Wigner crystal which is 1.1061. 
To test the accuracy of Eq. Q for different values of 7 we have simulated the distribution of 
charges on the surface of a sphere and calculated their electrostatic energy. To speed up the 
simulations, instead of using the full continuum algorithm for the surface Thomson problem, 
we have studied its discrete version. Thus, a number of sites were randomly placed on the 
surface of the sphere with a uniform probability density. The charges were then restricted 
to move only on these sites. There is a significant gain in the simulation time since all the 
electrostatic interactions can be tabulated once at the beginning of the simulation. 

In figure 1 we compare the result of numerical simulation with the analytic expression 
given by Eq. Q). The agreement is perfect for the whole range of n. 

Since the number of metastable states grows exponentially |10) with increase in n, our 
simulation is not able to locate the exact ground state for large number of charges. However, 
the energy of nearly degenerate metastable states is very close to that of the ground state, and 
the error thus accrued is minimal. This can be clearly seen from the absence of any visible 
fluctuations in the data points plotted in Fig. 1. 

We are now in position to compare the electrostatic energy of the configuration in which 
all the electrons are at the surface, with the configuration in which n — 1 particles are on the 
surface and one charge is at the center of the sphere, 

AEin) = E„_i + - En . (5) 
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Fig. 1 - Energy En for charges with 7 = 1.4, Eq. 0] as a function of the number n of particles, 
compared with the simulation data (points). The Madelung constant ~ 1.78 gives a perfect fit to 
the data points over the whole range on n. Were used 2000 sites for allowed positions of the charges 
on the surface of the sphere. 

From Fig. 2 we see that for 7 > 1, AE{n) starts positive, so that the charges are driven to 
the surface. However, as the surface particle population increases, lS.E{n) reaches a maximum 
and begins to dechnes. At the threshold number of charges n = nd"/) it becomes energetically 
favorable for the particles to penetrate into the bulk. For 7^1+ the critical number of charges 
diverges as. 

It is very curious that the Coulomb law is precisely at the border Hne of the two regimes. 
For the potentials with 7 < 1 repulsion is strong enough to drive all the charges to the 
surface, while for 7 > 1 and n > ndj) the charges will penetrate into the bulk. Unlike for 
Coulomb charges, for particles interacting by potential, the Thomson problem requires a 
full three dimensional analysis of the distribution of particles inside the sphere. The generalized 
Thomson problem is, therefore, much more complex than its classical counterpart. 
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Fig. 2 - The difference in energy AE{n) between the configurations in which one particle is at the 
center of a sphere and n — 1 particles are at the surface, and the configuration in which all n charges 
are on the surface. The values of 7 range from 1 for the topmost curve, to 1.4 for the bottommost 
curve, in intervals of 0.1. The topmost curve 7 = 1 is a monotonically increasing function of n, while 
for 7 > 1 the curves after reaching a maximum decline. 
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